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Abstract. We present a comprehensive numerical investigation of non- universal 
parameters and corrections related to interface fluctuations of models belonging to 
the Kardar-Parisi- Zhang (KPZ) universality class, in = 1 + 1, for both flat and 
curved geometries. We analyzed two classes of models. In the isotropic models the 
non-universal parameters are uniform along the surface, whereas in the anisotropic 
growth they vary. In the latter case, that produces curved surfaces, the statistics must 
be computed independently along fixed directions. The ansatz h — Vact+ (rt)^/'^x + ^i 
where x is a Tracy-Widom (geometry-dependent) distribution and is a time- 
independent correction, is probed. Our numerical analysis shows that the non- 
universal parameter F determined through the first cumulant leads to a very good 
accordance with the extended KPZ ansatz for all investigated models in contrast with 
the estimates of F obtained from higher order cumulants that indicate a violation of 
the generalized ansatz for some of the studied models. We associate the discrepancies 
to corrections of unknown nature, which hampers an accurate estimation of F at 
finite times. The discrepancies in F via different approaches are relatively small but 
sufficient to modify the scaling law i^^/'^ that characterize the finite-time corrections 
due to T]. Among the investigated models, we have revisited an off-lattice Eden model 
that supposedly disobeyed the shift in the mean scaling as and showed that there 
is a crossover to the expected regime. We have found model-dependent (non-universal) 
corrections for cumulants of order n > 2. All investigated models are consistent with 
a further term of order t^^/"^ in the KPZ ansatz. 



1. Introduction 

Dynamics of self-affine surfaces is a fascinating topic of non-equilibrium statistical 
physics where pattern formation, stochasticity, scale invariance and universality are 
joined in a unified framework P, |2]. The KPZ equation [3] 

f = -v^/^ + ^(v/.f + e, (1) 

proposed by Kardar, Parisi and Zhang in 1986, emerged as one of the most prominent 
theoretical systems and established a universality class which encompasses a plenty of 
models [H This equation describes the non-conservative evolution of an interface 
h{x,t) subjected to a white noise ^ defined by (^) = and {$,{x,t)^{x' ,t')) = DS{x — 
x')S{t — t'). The KPZ universality class goes beyond the surface dynamics including 
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theoretical studies in random polymers [1] and fluid transport [5]. The understanding 
of the KPZ universality class has undergone a noticeable progress during the last few 
years with its experimental realization P, [TJ El E] as well as the achievement of analytical 
solutions m [101 [m ini [131 [H] of the KPZ equation in 1+1 dimensions. 

Self-affine interface evolution is featured by a dynamical regime where the interface 



width w, deflned as the standard deviation of the surface heights, w = y {K^) — {h)^, 
increases in time as a power law w ^ f^, and a stationary regime where the interface 
width depends on the system size as w ~ L". /3 and a are the growth and roughness 
exponents [1], respectively, that assume the exact values (3 = 1/3 and a = 1/2 for 
KPZ class in (i = 1 + 1 [3J. The universality in KPZ systems includes other important 
universal quantities related to the height distributions (HDs) [71 |T5]. Motivated by 
analytical determination of the HDs of some models of the KPZ class, namely, the 
asymmetric exclusion process [I6] and polynuclear growth model p!7l [T8] , Prahofer and 
Spohn [T7] conjectured that the heights for any KPZ system are asymptotically given 
by 

h = Voot + sx{rty^''x (2) 

where s\ = sgn(A), foo and F are non-universal parameters while x is a stochastic 
variable with universal properties. Moreover, the conjecture states that, in the 
dynamical regime, x depends on the growth geometry and is given by Tracy- Widom 
(TW) distributions from random matrix theory [TP] for flat and curved growth 
geometries. The Gaussian orthogonal ensemble (GOE) is expected for the former 
while Gaussian unitary ensemble (GUE) is expected for the latter. Therefore, the 
KPZ universality class splits into subclasses depending on the growth geometry. This 
conjecture has been verifled by many analytical results [H [TOl [Til [131 [H]) computer 
simulations [201 1211 ESI [23] and, most excitingly, in experiments [HI [71 [8l|9]. 

The TW distributions are conjectured for the asymptotic regime, but many 
experimental [HI [S], analytical [TUl and numerical f2T[ [22] works have shown the 
existence of a slow convergence with apparently universal properties. More speciflcally, 
the rescaled distribution P{q), where 

h v^t . . 

has a shift in relation to the TW distributions that vanishes as 

(g) - (x) ~ t-'/\ (4) 

An exception to the correction was obtained in the partially asymmetric simple 
exclusion process (PASEP) for the asymmetry parameter = 0.7822787862... 
The PASEP result suggests that corrections different from would appear in very 
particular situations. However, Takeuchi [23] reported off-lattice simulations of a 
radial Eden model where the flrst cumulant of the scaled distributions converges to 
the asymptotic GUE value as in contrast with other Eden versions [21j. It is 

somehow intriguing that a stochastic model without control parameters does not have 
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the usual correction. Still in reference [21], a deterministic shift rj was found in the 
polynuclear growth model (PNG) and in the totally/partially asymmetric exclusion 
processes (TASEP/PASEP), generalizing equation g to 

h = v^t + sx{TtY'\ + r] + ... (5) 

The correction in the first cumulant derives directly from equation It was also 
shown that higher order moments have corrections 0{t'''^/^) [23] and, consequently, all 
higher order cumulants of order n>2 decays as or faster. Sasamoto and Spohn 

found out a result similar to equation ([s]) in their solution of the KPZ equation [TOl [12] 
with a difference that rj is random, but still independent of x- 

Many models belonging to the KPZ class are analytically intractable with our 
current knowledge and Monte Carlo simulations become the best accessible method to 
probe the universality of their interface fluctuations. The analysis of flat growth can 
be efficiently performed using lattice models [221 EH] whose computer implementations 
are commonly simple, relatively quick and low memory demanding. Eden model [26] 
and its variations [2] are basic examples of radial growth belonging to the KPZ class. 
However, it is well known that radial growth in lattices is distorted by anisotropy 
effects [23 [2H1 [2ni EO]. Therefore, off-lattice Eden simulations [211 [2Sl EI] or suitable 
tricks using tuning parameters [211 EO] are needed to investigate interface fluctuations in 
the entire surface. Although off-lattice simulations of Eden growth are quite affordable 
in 1 + 1 dimensions [32], the generalization to higher dimensions is cumbersome [33]. 

In this work, we present extensive simulations of models belonging to the KPZ 
universality class to address the validity of the extended KPZ conjecture including the 
correction rj given in equation (|5]). We investigated models exhibiting either isotropic 
or anisotropic growth. In the isotropic growth the velocity is the same for all parts of 
the interface such that the entire proflle can be used to perform statistics. On other 
hand, in an anisotropic growth, in which the interface velocity varies along the interface, 
we must analyze fixed directions independently and, consequently, a large ensemble of 
samples are required since one or a few points from each sample are used for statistics. 

We have found that the procedure to determine the model-dependent parameter F 
via second cumulant may be troublesome due to strong and/or puzzling corrections to 
the scaling. We adopted an alternative method using the first cumulant derivative, that 
exhibits a monotonic convergence to the asymptotic value with a correction in all 
investigated models. This method allowed to determine shifts in the first cumulants 
always consistent with a decay reinforcing the generality of the KPZ ansatz. We 

also investigated the convergence of higher order cumulants and identified complex and 
non-universal behaviors. Finally, an additional term t~^^^ in the KPZ ansatz was found 
in all models considered here. 

The paper is organized as follows. Section |2] presents the analysis of fiat growth 
models concomitantly with the description of the numerical recipes applied in the 
present work. Sections |3] and |4] are devoted to isotropic and anisotropic radial growth, 
respectively. Section [5] presents the analysis of the droplet growth. We draw our 
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2. Isotropic flat growth 



We consider two versions of the restricted solid-on-solid (RSOS) model [31] and the 
ballistic deposition (BD) model pQ on initially flat substrates. Averages were computed 
using 10^ samples and system sizes L = 2^^ resulting in more than 2 x 10^ data points 
for the statistics. 

The RSOS model is defined as a random deposition obeying the height restriction 
Ah = \hj — hj^i\ < m. In the BD model, particles move normally to the substrate and 
permanently attach to the first nearest neighbor contact with a previously deposited 
particle. In both models, periodic boundary conditions are assumed and the time is 
increased by At = 1/L for each deposition attempt. 

Accurate estimates of the non-universal parameters are imperative to a reliable 
characterization of the corrections. The asymptotic velocity is accurately determined 
taking the time derivative of the mean height in equation ([s]): 



{h)t 



+ 



(6) 



Then, the velocity Voo can be extracted by extrapolating {h)t versus in a linear 
regression for t — t- oo. The inset of figure [T] shows a typical plot used to determine v^o 
for RSOS model with m = 1. Interface velocities for all investigated isotropic models 
are shown in table [U 

The parameter F, in terms of the constants of the KPZ equation, is given by 
r = ^A^|A| with A = D/u |T5l [8]. The parameters A and A can be determined 
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Figure 1. Numerical determination of the parameter F for the RSOS model grown 
on flat substrates, with Ah < I, using different cumulants. Solid lines arc non- linear 
regressions. The inset shows the interface velocity against t^^/-^ (symbols) and the 
linear regression used to determine the asymptotic velocity. 
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model 




s\ 


Voo 


Ti 


r2 




r4 




RSOS (m= 




-1 


0.419030(3) 


0.252(1) 


0.2532(3) 


0.2525(2) 


0.2534(3) 


-0.32(4) 


RSOS (m= 


=2) 


-1 


0.60355(1) 


0.812(2) 


0.816(9) 


0.811(2) 


0.815(2) 


-0.66(6) 


BD 




1 


2.13983(1) 


4.94(1) 


4.778(2) 


4.799(7) 


4.74(2) 


-0.9(1) 


Eden D 




1 


0.51371(2) 


1.00(1) 


w 1 


0.99(1) 


0.98(2) 


0.50(5) 



Table 1. Non- universal parameters for distinct isotropic growth models. The number 
in parenthesis represents the uncertainties obtained from the non-linear regressions. 
The parameter r2 does not have uncertainty for off-lattice Eden model due to the lack 
of a monotonic convergence. 



numerically using the tilt dependence of the growth velocity and the amplitude of a two- 
point correlation function, respectively [35]. Alternatively, F can be obtained directly 
from equation ^ assuming that the cumulants of x ^'^^ those of TW distributions. The 
scaled second cumulant 

92=^^-^T^"{x\. (7) 

where denotes the nth cumulant of X, has been used to this purpose [HI ESI El [22] • 

In principle, one can obtain F from any higher order cumulant since equation ^ yields 

9n = f^,^r''l^r')c. n>2 (8) 
Alternatively, the parameter F can also be obtained from the first cumulant since 

g, = sy^?,m,-v^)ei'-^T'l\x)- (9) 

Although widely used in many experimental [HI EJ [7] and computer [221 EH] studies, 
the determination of F using second order cumulants has complications when there are 
statistical dependencies among x, r] or other unknown terms in equation ([s]). In this 
case, crossed terms appear in cumulants leading to relevant corrections that, in principle, 
may be puzzling. Otherwise, the analysis using {h)t is free from crossed terms and does 
not depend on rj. Noise in numerical derivatives counts against this last method. We 
propose that the first cumulant derivative yields more reliable estimates of F than {h?)c 
does, as discussed in the this and in the following sections. 

We denote by F„ = [gn/ (x")c)^^"'5 the value of F estimated via the nth cumulant. 
Figure [T] shows the curves used to determine the non-universal parameters of the RSOS 
model with m = 1. Independently of the used cumulant, the data have finite-time 
corrections. To estimate the asymptotic value as well as the scaling of the correction, 
numerical values of F„(t), that are recorded exponentially spaced in time, are plotted 
as function of s = Int. Assuming that there is a power law correction, we performed 
non-linear regressions in the form 

Tn = F„(oo) + cexp(-a„s), (10) 

discarding very short times. The extrapolated values for F„ are shown in table [T] 
All cumulants yield essentially the same asymptotic value for F„, that is in very good 
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agreement with our previous estimate using longer growth times (and less samples) 
without extrapolations [22]. The correction in Ti has an exponent ai = 0.661 ^2/3. 
Assuming that the next term in equation ([s]) is a power law C,t~'^ , we have that 

gi = T^"{x) + {Ot-^-"'. (11) 

So, a correction ai = 2/3 in (71, or equivalently in Fi, shows that the next leading term 
in equation ^ decays as t^^/^. The correction in that have the same exponent as 
{q^)c — {'}d)ci is 02 ~ 1 while 03 ^ 04 2/3 are found for the third and forth order 
cumulants. These results are corroborated in figure |2| where the differences between g^, 
n = 1, .., 4, and their asymptotic values (see table [2]) are plotted against time. It is worth 
mentioning that corrections in cumulants were already studied for RSOS model [22] and 
are consistent with the present results. 

The power law correction t~^^^ in the mean is also verified with a good precision as 
shown in the inset of figure [2| where the difference between the scaled height and the 
GOE first cumulant has a plateau when rescaled by the expected power law. Considering 
the variation of the plateau inside errors of F and foo, the estimated amplitude of the 
correction is (r^) ~ 0.32(4). The mean (77) for all isotropic models are shown in table [l] 

We also simulated the RSOS model allowing a height difference A/i < 2. The 
results are essentially the same except by a non-monotonic time-dependence for r2 that 
initially decays and then increases towards the asymptotic value also with a correction 
faster than For other gm the corrections are approximately t"^/'^, the same as in 

the m = 1 case. The non-universal parameters for m = 2 and cumulant ratios are also 
shown in tables [T] and [2l 

Figure [3] shows the curves F„ against time and the respective non-linear regressions 




Figure 2. Corrections in g„ for the RSOS model grown on fiat substrates. Daslied 
line represents a power law and the solid ones represent the power law t^"^/^ . Inset 
shows the difference between the first cumulant of equation ([3| and the GOE mean 
(x) = -0.76007 using Voo = 0.419030 and T = 0.252. The shift is scaled by the 
theoretic power law t^^/^. 



Non-universal parameters, corrections and universality in KPZ growth 



7 



used to obtain their asymptotic values for the BD modeL Surprisingly, the extrapolation 
from Fi is a little larger than the others. The difference of 3-4% is sufficient to affect 
the scaling law ruling the shift of the distribution in relation to GOE as shown in the 
inset of figure [3j The value F = 4.94 obtained from the first cumulant derivative yields 
an excellent agreement with a correction while the values of F extrapolated from 
higher order cumulant s do not. 

The cause of this difference is unknown. A possibility could be crossed terms 
involving random variables in higher order cumulants, that are not present in the first 
one. These terms would introduce complicated corrections leading to a non-monotonic 



convergence and F could not be properly extrapolated using equation (10) for n > 2. 



Another possible explanation is that the asymptotic distribution of q for BD model 
does not converge to GOE but to a shifted GOE distribution with x + where a 
is a deterministic and non-universal parameter. In this case, the deterministic shift 
a cannot be obtained using higher order cumulants since (a")c = for n > 2. If 
this last hypothesis is correct, we would have Fi = F(l + a/{x)Y that, according 
to the parameters of table [1} provides a tiny shift a ~ 0.01. Such a small shift is 
easily unnoticed in scaled height distributions obtained in simulations or experiments. 
Indeed, the correction t~^/^ shown in the inset of figure [s] is equally obtained in curves 
(9) ~ (x + 'j^) against t (data not shown). It important to emphasize that this shift 
was not predicted in the other KPZ systems investigated either in the present work or 
elsewhere. Independently of the numerical evidence, we believe that the non-monotonic 
convergence is the most likely scenario. 

In reference [22], we have showed that the shift in the first cumulant is consistent 
with the usual t~^/^ decay using the parameter F2 = 4.90(5) obtained from (/i^)c in 
simulations up to t = 5 x 10^ without extrapolation. We checked that the extrapolation 
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Figure 3. Determination of the non-universal parameter F for the BD model grown on 
flat substrates. Solid lines are non-linear regressions. Inset: Shift in the first cumulant 
scaled by the expected power law i^^/'^ using two estimates of F and Wqo = 2.13983. 
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of the data of reference [22j is consistent with our current estimate of A careful view 
in the double-logarithmic plot presented in reference [22] reveals a slightly bent curve 
indicating a correction not so close to t~^^^ as wished. Also, the uncertainties in F and 
Voo reported there are much larger than the current ones such that the propagated error 
in the power law is sufficiently large to embrace the exponent -1/3. 

Corrections in gn for BD model are shown in figure |4j In (71, we have found an 
exponent ai ~ 2/3, the same value found for the RSOS model. Nevertheless, an unusual 
exponent 02 ~ 1/2 was observed in second order cumulant. To our knowledge, there is 
not analytical or experimental analogues for this correction. This exponent is at odds 
with our previous analysis reporting a decay faster than in plots {(f)c—{'X^)c against 
time [22] . The discrepancy is due to the over-estimation of r2 used in the rescaled height 
q while in the current analysis, the dependence on r2 is implicit to g2 and, therefore, 
does not involve error propagation. The third and fourth order cumulants decay with 
exponents 03 ~ 0.9 and 04 ~ 1.25, respectively, that are consistent with the previous 
report 0{t~'^/^) or faster [22]. These results show that the exponents measured in the 
power laws related to the finite-time corrections are very sensitive to F and, thus, a 
precise characterization of the corrections demands accurate estimates of F. 




Figure 4. Corrections in g„ for the BD model grown on fiat substrates. Tlie lines 
are power laws i^^/s ^ f-i/2 ^ 2), {n = 3), and i^^ ^^ (n = 4). 



3. Isotropic radial growth 

We simulated radial growth using an Eden model proposed in reference [23] , that we have 
named as Eden D [33] since A, B, and C versions are already defined in literature [2|[30]. 
The model is defined as follows: In each time step, a particle of the cluster and a position 
tangent to it are randomly chosen. A new particle is added to the chosen position if this 
event does not imply the overlap with any other particle of the cluster. In the case of 
overlap, simulation proceeds to the next step. The time is increased by At = 1/A^, where 
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is the number of particles of the cluster at time t. Optimization strategies described 
in reference [32] were used to speed up the simulation. A cluster of diameter 8000 takes 
typically 8 min of simulation in a CPU Intel Xeon 3.2 GHz, whereas if no optimization 
is used the same simulation takes several hours of computation. We simulated the same 
version of the Eden model proposed in reference [23] using a much better statistics 
(90000 samples against 3000) and much longer growth times {t ~ 8000 against 2000). 
Considering the same interval used in reference [23] (250 < t < 2000), we have found 
Voo = 0.51390 in full agreement with Vqo = 0.5139(2) reported there. However, a more 
accurate estimate v^o = 0.51371(1) was found for t > 2000 (bottom inset of figure [s]). 

Curves used to determine r„(oo) for the Eden model are shown in figure [6] The 
second order cumulant depends non-monotonically on time as highlighted in the inset of 
figure|6} This non-monotonicity hampers an extrapolation to r2(oo) using equation (10). 
Using an extrapolation similar to that used in reference [23], i.e., forcing a correction 
02 = 2/3 in we obtained r2 ~ 0.995 considering the time interval 1000 < t < 8000. 
This value is slightly smaller than r2 = 1.02(2) reported in reference p3]. The first 
cumulant derivative yields a monotonic convergence with an exponent ai ^ 2/3 and 
extrapolates to Fi = 1.00(1). Moreover, and r4 also have corrections consistent with 
and converge to values very close to 1 (see table [l|. 

The main plot of figure |5] shows the shift (q) — {x) against time for the Eden model. 
Using the velocity Voo = 0.51390 obtained by Takeuchi, his result, a decay much faster 
than [22], is reproduced. If the more accurate estimate v^o = 0.51371 is used 

instead, a decay close to is observed for short times and is followed by a crossover 
to a smaller exponent. The top inset of figure |5] shows the local exponent against time. 
One can clearly see that the scaling exponent has not reached a stationary value even 




Figure 5. Correction in the first cumulant of the scaled height in relation to 
the corresponding GUE value (x) = —f. 77109 for off-lattice Eden D model. The 
parameters used were F = 1.00 and interface velocities indicated in the legends. Dashed 
lines are power laws t~^^^ and t~^^^ as guides to the eyes. Top inset shows the local 
exponent against time, while the bottom one shows the velocity against t^'^^^. 
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Figure 7. Corrections in scaled cumulants for Eden model. Dashed lines represent 
the power law t~'^l'^ as guides to eyes. 



for our longest simulated times. 

Corrections in are shown in figure [7} The quantity gi has an excellent agreement 
with a power law resulting again in a correction t~^^^ in the KPZ ansatz given 
by equation ([s]). Thus, a higher order term is expected in (g) — (x) and explains 
the short-time decay observed in figure [5j Performing a double power law regression 
{q) ~ {x) = at~^/^ + bt^"^^^, an excellent fit is obtained with positive amplitudes 
a = 0.50(5) and b = 4.0(1). This result shows that the amplitude of the correction 
rj is much smaller than the next leading term resulting a very slow crossover to 
the asymptotic scaling law t~^^^. Due to the lack of monotonicity, scaling of the 
corrections in g2 could not be determined with the present data. The correction in 
the third cumulant is Despite of large fluctuation, the forth cumulant also decays 
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consistently with 

Even though all parameters in table [T] depend on the non-universal parameter F, 
they can be combined in dimensionless cumulant ratios that must reflect the universality 
of X- The relative variance R, the skewness S and the kurtosis K defined by 



and 



R 



S 



K 



92 ^ 

9l~ 

93 
3/2 
92 



3/2' 

c 



9 A {X^)c 



9l 



ix') 



(12) 



(13) 



(14) 



must therefore be model-independent. 

Table [2] shows the cumulant ratios for all isotropic models investigated. The results 
are in excellent agreement with the corresponding GOE {Rgoe 



1.1046, Sgoe = 0.2935, 



and K, 



goe 



0.1652) for flat models and GUE [R. 



gue 



0.2592, S, 



gue 



0.2241, and 

Kgue = 0.09345) for the radial Eden model. Scaled height distributions P{q) for all 
isotropic growth models were reported elsewhere [221 [23] and, therefore are omitted 
here. We present scaled HDs for anisotropic growth models in the following sections. 



model 


91 


92 


53 


.94 


R 


S 


K 


RSOS (m=l) 


-0.4795(5) 


0.2553(1) 


0.03772(5) 


0.01075(5) 


1.110(2) 


0.292(1) 


0.1649(5) 


RSOS (m=2) 


-0.709(1) 


0.5547(4) 


0.1216(2) 


0.0512(1) 


1.103(4) 


0.294(2) 


0.1664(5) 


BD 


-1.294(1) 


1.8095(6) 


0.712(1) 


0.536(3) 


1.080(2) 


0.292(2) 


0.163(1) 


Eden 


-1.770(2) 


wO.81 


0.161(2) 


0.0607(5) 


0.26 


0.22 


0.093 



Table 2. Asymptotic values for gn for isotropic growth models. Dimensionless 
cumulant ratios R = 52 /.9i, S = 53/52^^ K = g^lg\ are also included. Lack 
of uncertainty in gi did not allow to determine errors in dimensionless ratios for Eden 
model. 



4. Anisotropic radial growth 

Radial interfaces were simulated using Eden models on square lattices. Two definitions 
are necessary to describe the models. Periphery sites are occupied sites with at least 
one empty nearest neighbor (NN), while growth sites are empty sites with at least one 
occupied NN. We investigated two versions. In the version Eden A, a growth site is 
chosen at random and occupied. For each choice, the time is increased by l/A^, where 
Ng is the number of growth sites. In the version Eden D on a lattice, a periphery site 
and one of its NNs are randomly chosen. If the selected NN is empty, it receives a new 
particle, otherwise, the simulation runs to the next step. For each attempt, the time is 
incremented by 1/Ap, where Ap is the number of periphery sites. In all versions, the 
growth starts with a single occupied site at the center of the lattice. 
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Figure 8. (a) Successive borders of on-lattice Eden D model for small sizes. The 
circle has the mean radius of the border, (b) First quarter of the border of a large 
on-lattice Eden D cluster. The inset shows a zoom of the cluster top. 

The start point to check the equivalence between isotropic and anisotropic growth 
is the interface width w, defined as the standard deviation of the radius measured in 
relation to the center of the lattice. At early times, it behaves as a power law w ~ t'^, 
where the exponent /3 = 1/3 is expected for the KPZ class m. d = 1-1-1. However, it 
is well know that Eden clusters are affected by the lattice-imposed anisotropy [271 EO] ■ 
More precisely, the axial direction (10)[|] grows slightly faster than the diagonal (11) 
implying that the dispersion around the mean radius of the border is asymptotically 
ruled by a diamond-like shape that produces a crossover from w ~ t^/^ at short times 
to w ~ t at long times [23 EO] • Figure [s] shows the evolution of an Eden cluster where 
the faster growth along (10) can be seen. 

Since different directions have different growth velocities, we must focus on the 
fiuctuations along a fixed direction. The radius on a given direction is defined as the 
distance from the origin of the farthest cluster particle along that direction. Equation (|2]) 
can be applied for an arbitrary direction with the parameters f oo and F depending on it. 
At a time t, we have a collection of radii r2{t), ■ ■ ■ rj^it)} along a given direction 

that are obtained from an ensemble of N independent simulations. We analyzed axial 
and diagonal directions that correspond to fastest and slowest growth directions in 
square lattices, respectively. We simulated 2.5 x 10^ clusters of diameter 5 x 10'^. 
Therefore, our statistics is performed with 10'' points for each analyzed direction. A 
typical simulation of Eden A and D takes about 3 s and 7 s, respectively, in a CPU Intel 
Xeon 3.20GHz. 

Figure [9] shows the interface width against time for Eden A computed for distinct 
directions. Interface width along direction (10) is larger than along direction (11) {Aw ~ 
6%), but they follow scaling laws in time with growth exponents /3io = 0.335(7) and 

= 0.332(7), respectively, in excellent agreement with the KPZ universality class. 
Growth exponents for Eden A and D are shown in table |4j The inset B of figure |9] 

I Here, we have borrowed the crystallographic notation where (10) represents the directions of unitary 
vectors x,-^, y and -y. 
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shows the analysis to determine the interface velocity for Eden A simulations. A higher 
velocity along axial direction is evident. 




Figure 9. Interface width against time for on-lattice Eden A model. Inset (A): 
Effective growth exponent obtained from the derivative of Inw versus Int. Inset (B): 
Interface velocity against t~'^/'^ . 
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Figure 10. Determination of the non-universal quantity r^/'^(x) for on-lattice Eden 
models. Solid lines are non-linear regressions used to extrapolate gi. The inset shows 
the correction in gi against time. Dashed line has slope -2/3. 



Corrections in gi are shown in figure 10 and are consistent with a decay 
implying that the next leading term in equation ^ is in analogy with isotropic 

growth models described in sections |2] and [3j The scaled cumulants against time for 



Eden A are shown in figure 11 Along the direction (10), we have found corrections 
consistent with t~^ for the second cumulant in both Eden A and D. Along direction 
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Figure 11. Scaled cumulants — (r'")c/(s^i"/'^) for on-lattice Eden A afong 
directfons (10) (open symbols) and (11) (filled symbols). Lines are extrapolations 
to i — > oo. 



(11), g2 has a slight non-monotonicity that comphcates the extrapolation to long times. 
It was not possible to accurately resolve the scaling of the corrections for higher order 
cumulants but we have found that they do not decay slower than as illustrated 

in figure 12 This indicates that the correction rj is statistically independent of x since, 
otherwise, corrections would decay as for all higher order cumulants [B1I21]- Notice 
that many features of our current on-lattice simulations were observed for the isotropic 
growth, where was found for gi, g^ and g^ and a non-monotonicity for g2. Table [s] 
shows the non-universal parameters Voo and gn obtained in our simulations. 

One can see in table |3] that the estimates for Fi are slightly larger than those for 
r2 for all investigated models, as observed for ballistic deposition on flat substrates. In 
the rest of this work, we use gi to investigate the scaling and universality of the models 
since it yields reliable estimates of F as discussed in section |2} 



Table 3. Non-universal parameters of several anisotropic growth models with curved 
geometry. Acronyms: EdA (Eden A); EdD (Eden D); WRSOS (wedge restricted solid- 
on-solid); BDD (ballistic deposition droplet). The missing uncertainties for EdA/(ll) 
and EdD/(ll) are due to the lack of monotonicity in the second cumulant of the 
respective data. 



Model 


Voo 


9i 


92 


.93 


54 


Ti 


r2 


iv) 


EdA/(10) 


2.1824(3) 


-2.927(3) 


2.195(1) 


0.728(2) 


0.45(1) 


4.51(1) 


4.43(1) 


1.7(1) 


EdA/(ll) 


2.1401(3) 


-2.807(3) 


1.956 


0.598(2) 


0.345(9) 


3.98(1) 


3.73 


1.1(1) 


EdD/(10) 


0.6186(1) 


-1.754(2) 


0.779(1) 


0.1535(8) 


0.057(1) 


0.971(3) 


0.938(2) 


1.6(1) 


EdD/(ll) 


0.61055(5) 


-1.677(2) 


0.695 


0.127(2) 


0.044(1) 


0.85(1) 


0.79 


1.1(1) 


WRSOS 


0.41903(1) 


-1.115(4) 


0.3233(2) 


0.0418(2) 


0.0095(4) 


0.250(3) 


0.251(1) 


-0.7(1) 


BDD 


2.13984(2) 


-3.023(9) 


2.331(5) 


0.802(6) 


0.49(1) 


4.97(4) 


4.85(2) 


1.3(1) 



Due to the anisotropy, the surfaces studied here do not have translational invariance 
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Figure 12. Corrections in the second cumulant for on-lattice Eden models along 
direction (10). The dashed line is a power law with exponent —1. Inset: Corrections 
in the third cumulant. The dashed lines are power laws with exponents —2/3 and — 1. 



Table 4. Universal quantities for several anisotropic models with curved geometry. 
Acronyms as in table [sj Cumulant ratios are defined by equations ( 12 )-( 14 1. The KPZ 



cumulant ratios correspond to the CUE distribution were taken from reference [17]. 
The lack of uncertainties in the second cumulant of EdA/(ll) and EdD/(ll) models 
did not allow the error propagation in the cumulant ratios. 



Model 


/3 


R 


S 


K 


EdA/(10) 


0.335(7) 


0.256(1) 


0.224(1) 


0.093(2) 


EdA/(ll) 


0.332(7) 


0.25 


0.22 


0.090 


EdD/(10) 


0.338(9) 


0.253(1) 


0.223(2) 


0.094(2) 


EdD/(ll) 


0.334(10) 


0.25 


0.22 


0.091 


WRSOS 


0.3326(5) 


0.260(2) 


0.227(1) 


0.091(4) 


BDD 


0.335(2) 


0.255(2) 


0.225(2) 


0.090(2) 


KPZ 


1/3 


0.2592 


0.2241 


0.09345 



so that two-point local quantities are not well-defined. So, the parameter A and, 
consequently, F cannot be directly measured if cumulants of x ^^^e unknown. We then 
define 

, r-Voot _ q , . 

= 

such that (g') — i- 1 for t — )■ oo. We can obtain the scaling law ruling the shift in the 
mean by plotting 1 — (g') against time using only the directly measurable parameters 
Voo and gi. The results are shown in figure 13, where we have obtained corrections very 
consistent with already observed in many other systems |S1 El HDl Ell [M] ■ The slow 
crossover to t"^/^ observed for off-lattice simulations of Eden D [23] was not observed 
in the on-lattice version. Moreover, our previous off-lattice simulations of Eden B[§] also 

§ The difference between Eden D and Eden B is that any nearest neighbor can be chosen in the former 
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yielded a shift consistent with t~^/^ [21j. An important remark is that the directly 
measurable amplitude in the scaling 1 — (g') ^ —{v)/9i^~^^^ ^ilso yields a way to obtain 
an estimate of {rj). The values of {rj) for anisotropic growth models are shown in table [s] 



o EdMlO) 

• EdMll) 
□ EdD/<10> 

* EdD/<ll> 
o WRSOS 
A BDD 




10- V 
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10' 



t/t 



10-' 
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Figure 13. Scaling law describing the shift in the mean of scaled surface distributions 
for several anisotropic KPZ models. The time was scaled by the maximum growth 
time reached in each model in order to improve visibility. The solid line represents the 
power law t"^/"^. Acronyms as in table [s] 

The dimensionless cumulant ratios obtained for both Eden models are shown in 
table 111 All ratios are universal and show remarkable agreement with GUE [T7] . Notice 



that using dimensionless cumulant ratios as equations (12), (13) and (14), all cumulants 



of X can be determined as functions of a single cumulant. This reasoning also applies 
to systems where x is unknown raising a general strategy to probe universality in the 
interface distributions. 



The top panel of figure 14 compares the distributions scaled accordingly 



equation (15) for Eden A in the direction (10) at different times with the GUE 
distribution scaled to a mean 1, x* = x/{x)- The distributions converge to GUE 
as time evolves but the shift is evident. If the estimated shift (77) = 1.7(1) is explicitly 
included by defining 

, r-Vopt- {r]) 

an excellent collapse between P{q*) and P{x*) is found (figure [l4| bottom). Similar 
results are found for all investigated models. 



5. Droplet growth 

A paradigmatic model for the growth of a droplet morphology is the PNG model, where 
islands are randomly nucleated over existing ones and, once nucleated, they starts 

while only the empty ones are eligible in the latter. 
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Figure 14. Top: Distributions for Eden A model at different growtli times scaled 
accordingly equation (15 1. Bottom: Same distributions rescaled using equation (16) 
that includes the shift (77) = 1.7. Insets show zooms around the peaks. 



to grow laterally with constant velocity. If we consider an initial island that grows 
indefinitely from a single nucleation at the origin, the PNG model has an exact solution 
with an asymptotic droplet shape, where the height at the origin is given by equation ^ 
with exactly known parameters [TTIIIH]. PNG model simulations with droplet geometry 
have been performed to be compared with analytic results [171 [20]. Other models and 
initial conditions leading to droplet interfaces has also been investigated [TOl [161 [36] . 

Droplet geometry can also be investigated in more complex models for which 
analytical results are currently unavailable as, for example, the BD model p8j. The 
growth of ballistic deposition droplets (BDD) starts with a single particle stuck to the 
origin of the system. Particles move ballistically along direction —y. If a particle visits a 
NN site of the aggregate, it irreversibly attaches to this position and becomes part of the 
aggregate. Let L{t) be the size of the active region defined as the set of x-coordinates 
where a growth can be tried. When a particle deposition is tried the time is increased 
by At = 1/L{t). Figure 15 shows a typical ballistic deposition droplet. 

The non-universal parameters depend on the position in the active region in analogy 
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Figure 15. Ballistic deposition droplet after deposition of 56000 particles. The first 
half of deposited particles are depicted in black and the rest in green. The solid line 
represents the interface. 



to the direction dependence in the Eden growth, whereas cumulant ratios and growth 
exponent do not. An ensemble of 2 x 10^ samples were used to compute statistics. 
Growth of a cluster for t = 10^ (~ 10^ particles) takes about 2 s in an CPU Intel Xeon 
3.20GHz. For sake of conciseness, we present only results for the height fluctuations at 
the origin. 

The curves used to determine the non-universal parameters (table Isl) are shown in 



figures 16 and 17 The interface velocity has the usual behavior illustrated in the inset 
(B) of figure 16 However, BDD growth has strong corrections in the growth exponent 
as confirmed in the inset (A) of figure 16, in analogy to its fiat counterpart [221 [37] . The 
correction in the mean converges in the usual way as shown in figure 



13 



Notice 

the presence of a crossover from a faster initial decay to the regime analogous 
to the crossover observed for the isotropic Eden D model. Similarly to the isotropic 
growth, corrections in g2 decays with an unusually small exponent t""'^^. 

The quantity gi approaches its asymptotic value as the same correction 

obtained for fiat simulations of isotropic growth models and for all versions of the 
anisotropic Eden model presented in section |4j The cumulant of order n = 3 also has a 
correction very close to t~^/^ while the cumulant of order ra = 4 has a correction faster 
than In particular, the correction observed in g^ is slower than that observed 

for the fiat case, demonstrating the intricate and non-universal scenarios involving 
corrections in cumulants of order n > 2. The inset of figure 17 shows the corresponding 
corrections in gn- Finally, the growth exponent and dimensionless cumulant ratios 
converge to the values expected for the KPZ class as shown in table |4j 

We also considered the restricted solid-on-solid (RSOS) growth model [31] with a 
wedge initial condition, the WRSOS model. At each time step, a site in the growing 
zone (of size L{t)) is randomly selected and its height increased by 1 if the constraint 
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Figure 16. Interface width against time determined at center of the BDD model. 
Solid line is a power law i^/'^. The inset (A) shows the effective growth exponent while 
inset (B) shows the interface velocity against t^'^/'^. 




Figure 17. Semi-logarithmic plot of gn against time for BDD model. Solid lines are 
the non-linear regression used to extrapolate the asymptotic values. Inset: Difference 
gn{t) — gnioo) against time. The dashed line represents a decay t~^^^. 



Ah = \h{j,t) — h{j ± l,t)\ < 1 is satisfied, otherwise, the deposition attempt is refused. 
The time is incremented by At = 1/L{t) for each attempt. A wedge initial condition 
h{j,0) = \j\, obeying Ah = 1, was considered. This initial condition imphes a droplet 
growth since the radius of the growing zone increases with average velocity 1 while in 
the center ((10) direction) it has a smaller velocity (Tab. [s]). Typical interfaces for 



distinct times are shown in figure 18 



Figure 19 shows the interface velocity at j = 0, 10, and 1000. They converge to 
the same asymptotic value observed for the flat geometry [22j , but differ at short times. 
This means a macroscopic shape, around the droplet center, that asymptotically moves 
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with constant velocity. The macroscopic shape is almost perfectly fitted by a parabola 
h{x) = h{0) + ax"^ in analogy with the solution of the KPZ equation with a sharp 
wedge initial condition [lOj. The surface statistics was computed using an ensemble 
of 5 X 10^ samples. A typical run up to t = 3000 takes about 0.2 s in a CPU Intel 
Xeon 3.20GHz. Non-universal parameters obtained for the height fluctuations at the 
origin j = are shown in table [3] The universal quantities shown in table |4] exhibit an 
excellent agreement with the KPZ universality class. Far from the origin, the accordance 
with the KPZ conjecture is also observed, but the farther from the origin the slower the 
convergence. 

The finite time corrections in the cumulants of WRSOS were also studied. The 



shift in the mean approaches the GUE one with the correction t (figure 13). The 



quantities for n = 1, 3 and 4 have a correction t while g2 decays slightly faster 



than t (figure 20) in analogy with the flat simulations of RSOS. 



1000 k 



500 - 




1000 



Figure 18. Typical interfaces obtained for WRSOS model for deposition times t = 512 
(lower) and t — 1024 (upper) . 



6. Discussion 



We have performed extensive simulations of models belonging to KPZ universality class 
in order to probe the generality of the KPZ ansatz given by equation ^ and to determine 
further corrections to this equation. We have analyzed two classes of models. In the 
isotropic growth models, all points of the interface are used to perform averages, in 
contrast with the anisotropic ones, where the non-universal parameters vary along the 
surface, limiting the statistics to a single or a few points per sample. 

During the numerical investigation, we have found out that the determination of 
the non-universal parameter P using the asymptotic value of the scaled second cumulant 
92 = {h'^)c/t^^'^ 1 has limited accuracy for some models, probably due to unknown 
puzzling corrections. We propose an alternative way to obtain P using the first cumulant 
derivative by means of the quantity gi = Sx3{{h)t — foo)^^^^ r^^^(x)- It is important 
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Figure 19. Interface velocity at different positions of WRSOS growth. The inset 
shows a zoom to show the asymptotic velocity. 
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Figure 20. Time corrections in gn for WRSOS model at j — 0. Dashed line is a 
power law t^"^/^. 



to mention that noise in the numerical derivative counts against this method. Here, 
we used very large statistics to smooth the curves. However, if statistics is limited, 
as possibly in case of experiments, one can use numerical procedures to smooth the 
derivative [39]. We also considered scaled cumulants Qn = {h'^)c/ {s'^t'^/^) of third and 
fourth order to determine F. In general, we have observed that small deviations in 
r may suggest a fake violation of the generalized KPZ ansatz. Estimates taken from 
gi give very good agreement with the generalized KPZ conjecture for all investigated 
models. Dimensionless universal cumulant ratios are not very sensitive to corrections 
in cumulants as one can see in tables |2] and [4], in which excellent agreements with GOE 
and GUE distributions were found for fiat and curved models, respectively. 
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Comparing the non-universal parameters for the fiat RSOS and BD models (table 
[T]) with its droplet counterparts (table [s]), we see that within the error bars they are 
equal, with exception of the shift {if). On the other hand, in the Eden D model all 
parameters are different for on- and off-lattice simulations. This is expected, since the 
lattice constraint should change the strength of fluctuations. 

In off-lattice simulations of the Eden model, we determined the shift in the height 
distributions and our results support a crossover from t^'^/^ to t^^^^ in the first cumulant 
in relation to GUE, clarifying an apparent exception to the correction t^^/^ recently 
suggested by Takeuchi |23] . 

From the numerical side, we propose that the estimates of F taken from the first 
cumulant are more reliable than those obtained from the higher order ones since, in 
principle, we know little about the nature of terms beyond x in the KPZ ansatz. These 
terms may introduce anomalous corrections as, for example, in the case of statistical 
dependence among random variables. In fact, in contrast to the correction in gi 
obtained for all investigated models, we have not found a standard in the corrections of 
higher order cumulants. 

The first cumulant derivative analysis for all investigated models are in agreement 
with an additional term generalizing the KPZ ansatz to 

h = v^t + s^iVtY'^ + ri + Ct~^/^ + ■■■■ (17) 

Corrections 0{t~^/^) were reported also for some analytically tractable KPZ class 
models [21], suggesting that this term is a general property of KPZ systems. However, 
in the solution of the KPZ equation with a sharp wedge initial condition a correction 
0{t^^) was reported [T2] . possibly because the non-universal parameter C, is zero in 
this case. 

The generalized ansatz given by equation (17) allows one to speculate on the origin 
of the complicated corrections obtained for cumulants of order n > 2. For instance, if 
rj and C, are statistically dependent, a term t~^ will appear in the second cumulant in 
addition to the leading term If the amplitude of t~^ is much larger than that of 

t ' , one would observe a slow crossover tot-2/3. This scenario explains the correction 
in g2 observed for the RSOS model. Still in the field of speculation, if x ^"^^ V ^^^^ 
statistically dependent the two first leading terms of the correction in g2 would be t~^^^ 
and A gain, if the amplitude of the second leading term is much larger than the 

first one, we could explain the anomalously slow exponent observed in the BD model as 
a crossover between the two scaling regimes. However, a double power law regression, 
as that used to determine the shift correction in Eden D, does not support such a 
conjecture. 

Now we discuss how universality in height distributions can be checked without 
an explicit knowledge of the parameter F. Performing a rescaling using the directly 
measurable parameters from tables [1} |2]and[3} we can define the variable: 

* ^ h-v^t- (rj) 



Non-universal parameters, corrections and universality in KPZ growth 



23 



The scaled distributions for this variable, P{q*), agree remarkably well with the GUE 
distribution scaled to unitary mean, x* = xl (x); shown in Figure 21 , top. Our entire 
analysis depends implicitly on (x), but does not provide information about it beyond 
its negative sign, since g\ < 0. We can alternatively define a variable with asymptotic 
unitary variance using the measurable parameter g2 rather than gi to find 

h - Voot - {r]) 



1 



X 



with a mean depending on the dimensionless and measurable cumulant ratio R given 
by 

(g**)=sgn((x))/?-'/' = -1.9641. 



Again, as shown in figure [21] bottom, a very good collapse is obtained. We remark 
that this kind of analysis could be very profitable to probe the universality of height 
distributions in systems where they are not known a priori^ such as high dimensional 
KPZ growth or models in universality classes without analytical counterparts. Although 
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Figure 21. Rescaled distributions for all investigated models using equations (16) 
and ( 18 ) in top and bottom panels, respectively. 
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the complete distribution are not determined, once (x) or y (x^)c is fixed everything 
concerning the HD can be derived. 
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